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The QCD phase transition is studied on 16'^ and 32'' x 4 lattices both with and without quark loops. We 
introduce a new zero-flavor or quenched species of quark and study the resulting chiral condensate, {^Q as a 
function of the ^ mass, m^. By examining {^(^) for 10"^" < < 10 we gain considerable information about 
the spectrum of Dirac eigenvalues. A comparison of ma = 0.01 and 0.025 shows little dependence of the Dirac 
spectrum on such a light, dynamical quark mass, after an overall shift in /3 is removed. The presence of sufficient 
small eigenvalues to support anomalous chiral symmetry breaking in the high temperature phase is examined 
quantitatively. In an effort to enhance these small eigenvalues, {^Q is also examined in the pure gauge theory in 
the region of the deconfinement transition with unexpected results. Above the critical temperature, the three Z3 
phases show dramatically different chiral behavior. Surprisingly, the real phase shows chiral symmetry, suggesting 
that a system with one flavor of staggered fermion at At = 4 will possess a chiral a phase transition — behavior 
not expected in the continuum limit. 



1. INTRODUCTION 

The spectrum of the Dirac operator is closely 
related to a number of important aspects of fi- 
nite temperature QCD. First, a non-zero den- 
sity of small eigenvalues is necessary for spon- 
taneous chiral symmetry breakingjl). Second, if 
the QCD phase transition is second order, one 
expects this eigenvalue spectrum to show critical 
behavior for T « Tc- Third, the Dirac spectrum 
should be particularly sensitive to the effects of 
quark loops since the Dirac determinant for small 
quark mass strongly suppresses gauge configura- 
tions with small Dirac eigenvalues. Finally, the 
universality arguments which predict the order of 
the QCD phase transition!^, ^'sly on the explicit 
breaking of the continuum Ua{^) symmetry by 
the axial anomaly. Such anomalous symmetry 
breaking requires a significant density of small 
eigenvalues. 

Since it is difficult to obtain the information 
about the spectrum directly we choose to exam- 
ine a less familiar quantity: the quenched or zero- 
flavor chiral condensate for an auxiliary quark 
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field (^{x) with mass m^. By studying the depen- 
dence on m,^, we obtain interesting, qualitative 
information about the spectrum, which suggests 
unexpected behavior near the QCD phase transi- 
tion. 

We begin with a normal, full QCD simulation 
at finite temperature described by the partition 
function: 

Z,n,p =UiJ dme-f"^^ det{D + m)^f/\ (1) 

Here the integral is performed over all link matri- 
ces Ui and the gauge action Sg is given by 

Sg = -\Y.retrUv (2) 

with U-p the ordered product of the four matrices 
associated with the fundamental plaquette V. We 
use the staggered fermion Dirac operator D. 

The determinant entering the partition func- 
tion in Eq. |l| represents the effects of quark loops, 
for Nf flavors of quarks of mass m. We will refer 
to the corresponding lattice quark fields as x{^)- 
In this paper we introduce a second species of 
fermions, C(^)j with mass tuq which will appear 
in observables but does not enter the quark de- 
terminant. These C fermions may be thought of 
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as quenched or TV/ = quarks. We then com- 
pute the expectation value of the quantity in 
the thermal background specified by the partition 
function of Eq. |l|. Thus, 



' {D + mc) 



\x) 



(3) 



where the average ( )^ „ is calculated using the 
distribution of gauge fields in the QCD path inte- 
gral. Although this quantity is no easier to com- 
pute than the normal chiral condensate (xx) j f he 
mass can be varied without changing the mass 
m which enters the Boltzmann factor and which 
must be fixed during a particular Monte Carlo 
run. Typically we study 25 values of m,^. 

This quantity is of special interest because it is 
a simple transform of the Dirac eigenvalue den- 
sity. If we evaluate the propagator in Eq. ^ by 
inserting a complete set of eigenfunctions of D, 
and averaging over x, we find 
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Here we have exploited the 7^ symmetry which 
makes the Dirac spectrum symmetric about zero 
and introduced p(A), the Dirac eigenvalue density 
per unit volume, averaged over the thermal dis- 
tribution of Eq. 1^. The function p{X) depends on 
f3 and the dynamical quark mass m. 

Thus, (CC) is a transform of the Dirac spec- 
trum, strongly weighted toward eigenvalues A ^ 
rriQ. Note the C/(l) chiral symmetry for the quark 
field C, present in the staggered fermion formula- 
tion, implies that (C^) should vanish as niQ —> 0. 
Spontaneous breaking of this symmetry results if 
p(0) 7^ 0. In which case the explicit factor of 
in the numerator of Eq. |^ is compensated by the 
linear divergence of the A integral as niQ 0. 
The result is the Banks-Casher formula||l[]: 



(CC)(mc = 0)=MO). 



(5) 



Thus, we argue that the function (CC)("^c) i^ 
easy to compute, contains information about the 
spectrum of Dirac eigenvalues and is closely con- 
nected to spontaneous symmetry breaking. In 
fact, as we discuss below, the infrared behavior 
of (CC)('^C ^ 0) also determines the character of 



anomalous chiral symmetry breaking, a quantity 
difficult to study directly with lattice fermions. 
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2 CHIRAL TRANSITION 



2.1. Critical Behavior 

Figure ^ shows the quantity {^0 as a function 
of rriQ for a series of values of f3 in the critical 
region. These results were obtained on a 16'^ x 4 
lattice with two degenerate flavors of dynamical 
fermions of mass ma = 0.01. Note, this is a log- 
log plot showing six orders of magnitude variation 
in m^. As can be seen, the behavior is quite com- 
plex. 
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Figure 1. The expectation value (CC) plotted as 
a function of the mass for a series of values of 
f3. Shown arc two-flavor calculations with quark 
mass ma — 0.01 on a 16"^ x 4 volume. 



First consider /3 in the low temperature, chi- 
rally asymmetric region, f3 < f3c- For m(^ > 
10^^, (CC) shows a complex non-linear behav- 
ior much of which is certainly due to finite lat- 
tice spacing effects. (In the continuum limit a 
simple linear behavior is expected for large 
with a quadratically divergent coefficient.) For 
10^^ < TO^ < 10~^ a constant plus linear fit 
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works well, at least for those values of (3 farthest 
from (3c- For < 10^'' finite volume effects 
set in and {C,Q) approaches zero linearly as ex- 
pected. It is interesting to note that a somewhat 
more elaborate fit, writing as the sum of a 
constant and fractional power, (^C) = co + ci 
works well over a significantly larger mass range 
10^'* < < 0.05 using a single value of y « 0.6. 
Such a fit merges nicely with a description which 
holds at /9 = /?c ~ 5.275, where only a frac- 
tional power, close to 0.6, is seen. Figure || 
shows our (3 = 5.245, 5.2625 and 5.265 data to- 
gether with these constant plus linear and con- 
stant plus fractional power fits. For the smallest 
value, (3 = 5.245, the simple constant plus linear 
fit works nicely. However, for the more nearly 
critical (3 = 5.265, Figure ^ shows the constant 
plus fractional power working visibly better. 

For (3 > 5.275, a fractional power, = cmfj., 
begins to fit well with the power growing with 
increasing (3 so that by (3 = 5.34 the power is 
y = 0.98, very close to 1.0, the expected linear 
dependence on at high temperature. Figure |^ 
shows these fits for f3 = 5.27, 5.275, 5.285 and 
5.34. For the smallest (3 value, 5.27, a pure frac- 
tional power does not fit well — a constant term is 
also needed while for the larger values of /? a pure 
fractional power, with zero intercept, works quite 
well. At these larger values of P the fits work at 
even lower masses suggesting that the region of 
small eigenvalues, A <« 10^"*, which may be dis- 
torted by finite volume effects, is less important. 

If we attempt to study the critical behavior of 
ixx) a function of the quark mass m, there are 
many predictions of universality For example, 
at Tc, one expects (xx) = '^i tt^^ for m ^ 0. 
Of course, this singular behavior in the mass can 
arise from a combination of the dependence of 
(xx) on the quark mass as it enters both the prop- 
agator and fermion determinant. This dual role 
played by the quark mass can be seen explicitly 
by relating (xx) and (CC) : 

{Xx)\l3,m ^ {COi'mc =m)\i3,m (6) 

Here the mass in the quark propagator explic- 
itly entering the expectation value, to<^, is often 
referred to as the valence mass, TOvai, while the 
mass in the fermion determinant, m, used when 
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Figure 2. An expanded version of the previous 
Nf — 2 plot showing possible fits for (3 < (3c- 
These (d values all show a non-zero intercept as 
niQ —>■ 0, when extrapolated from the region 
fnc > 5 X 10"**, a region apparently unaffected 
by finite volume distortions. Note how the con- 
stant plus fractional power fit agrees with the sim- 
ulation over nearly three orders of magnitude for 
P = 5.265. 



computing {^C.) is the sea quark mass, msea- 

Thus, it is the dependence on mvai that we are 
investigating when we examine (C^) as a function 
of . Only if the simulations do not depend in 
a critical way on the quark mass in the determi- 
nant, should we attempt to compare our exponent 
y\f3=f3, with 1/S- 

2.2. The High Temperature Phase 

The fact that (C^) appears to vanish with van- 
ishing m(; when /3 > /3c is a bit surprising. 
The conventional picture assumes the presence of 
topological configurations that break the ?7a(1) 
symmetry anomalously. This requires a signif- 
icant density of small eigenvalues in the Dirac 
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Figure 3. An expanded version of the previous 
Nf = 2 plot showing possible fits for (3 > (ic- 
These (3 values, except perhaps for /3 = 5.27, all 
show a zero intercept as — > 0, when extrapo- 
lated from the region mc^ > 5 x 10^''. Although 
a pure fractional power does not fit very well at 
the smallest /3 value shown, 5.27, for larger (3 this 
simple description works very well. 

spectrum. For example, a dilute instanton gas, 
expected by some to describe the region above 
the critical point, would give a non-zero density 
of eigenvalues, p(0) ^ 0. Recall, that it is only 
for the physical quantity (xx) that one expects 
the restoration of chiral symmetry. In the limit 
m — > not only is the valence mass in the prop- 
agator, TTivai vanishing, but also the mass in the 
determinant, TOgea is going to zero. It is the added 
suppression of small eigenvalues, coming from the 
zero quark mass limit of the fermion determinant, 
that is presumably required for chiral symmetry 
restoration above T^. 

Thus, a small sea quark mass should suppress 
the small eigenvalues associated with anomalous 
breaking of Ua{^) symmetry in the high temper- 
ature phase. Since in this calculation we do not 
set rrisea to zero, one might expect that anoma- 
lous breaking of symmetry will require 
{CQiraQ ^ 0) 7^ 0, contrary to what we see. Of 
course, this non-zero value of {(,(,) will be sup- 



pressed by factors of mseaO = 0.01. For that 
reason we will examine larger dynamical quark 
masses and even a quenched calculation below. 

2.3. Anomalous Ua{^) Symmetry Breaking 

Given the suppression of small eigenvalues that 
we see for /3 > /3c it is natural to ask if we can see 
any evidence for anomalous symmetry breaking. 
This is an important question given the effects 
of such symmetry breaking on the order of the 
two-flavor phase transition. 

Consider a continuum, two- flavor theory with 
SU{2) X SU{2) flavor symmetry (in the limit of 
vanishing fermion masses) above the phase transi- 
tion. The axial anomaly is expected to break the 
anomalous Ua{^), forcing the two-point functions 
constructed from the operators ijj^^T^ip and ipT'"ijj 
to be unequal. (Here ip is the continuum fermion 
field.) If these two-point functions are evaluated 
at zero momentum and subtracted, their anoma- 
lous 1/^(1) symmetry breaking difference (here 
called Lo) can be written as follows: 
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p{X)dX 

(A2 -I- to2)2 



(7) 



We can compare this equation with our results 
for the behavior of p(A) as A —> at fixed quark 
mass. In fact, the behavior of {(0 as for 

/3 just above Pc suggests that p(A) ^ c[m)\y for 
< y < 1 which would suggest lo ~ c{m) m^y~^\ 
However, the fermion determinant may well sup- 
press the small eigenvalues, with c(m) vanishing 
as m goes to zero. 

We can relate w directly to our quantity 
(CC)(mc): 
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dniQ 



(8) 



and attempt to study how it varies as m ^ 0. 
Figure ^ shows the quantity oj for a series of val- 
ues oi (3 > (3c at ma = 0.01 and at /3 = 5.327 for 
ma = 0.025. Clearly oj is non-zero in this region 
and does not appear to decrease as ma decreases 
from 0.025 to 0.01 provided we adjust (3 to remain 
a fixed displacement above (3c- Thus, this calcu- 
lation appears to show the presence of anomalous 
symmetry breaking for Nf = 2, in the high tem- 
perature phase, just above the transition. 
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Figure 4. The quantity lu which directly mea- 
sures the strength of anomalous symmetry break- 
ing plotted versus /3. On the left side results for 
ma = 0.01 are shown while on the right, results 
for ma = 0.01 and 0.025 are compared at /3 = 5.3 
and 5.327 respectively. 



/3c found in our 0.01 and 0.025 simulations. 



<^^> 0.25 



- 5.245 

- 5.255 ma=0.01 

- 5.26 ma=0.01 

- 5.272 ma=0.025 



i 

31 m 



3; S 



3. EFFECTS OF QUARK LOOPS 

In order to address the effects of the fermion 
determinant, let compare our ma — 0.01 and 
0.025 results in greater detail. We find that 
even when comparing the quantity {^C) evaluated 
over the entire range of m<^, the ma = 0.01 and 
ma = 0.025 simulations can be made to agree 
within a few percent if we allow for a quark-mass 
dependent shift in f3. In Figure || we show (CC) 
computed at /3 = 5.272, ma = 0.025 compared to 
a group of ma = 0.01 results. The quite precise 
5-6% agreement with the (3 — 5.255, ma — 0.01 
curve over the entire TO(^ range is striking. Sim- 
ilarly, for higher temperatures. Figure ^ shows 
(CC) computed a.t (3 — 5.327, ma ~ 0.025 com- 
pared to ma = 0.01 results. Again, we see 5-6% 
agreement with the /3 = 5.3, ma — 0.01 curve over 
the entire range. Thus, in both the chirally 
symmetric and asymmetric phases, the change in 
dynamical quark mass from 0.01 to 0.025 can be 
quite accurately compensated by a simple shift in 
p. No disparate effects on low or high eigenval- 
ues are seen. The shift in /3 of 0.027 needed in 
the high temperature phase is precisely the shift 
in (3c that we identified earlier when comparing 



Figure 5. The quantity (CC) plotted as a func- 
tion of m(; for p = 5.272 and ma = 0.025. It is 
compared with ma = 0.01 results for three values 
of f3. The close agreement with the ma = 0.01, 
P = 5.255 curve is evident. 



It is interesting to note that the 0.017 shift in 
P found in the low temperature regime is sub- 
stantially smaller than the 0.027 shift needed at 
high temperature. Thus, it is clearly incorrect to 
describe the effects of quark loops as causing a 
simple shift in p. At the least, that shift is itself 
/3-dependent. 

4. QUENCHED CHIRAL TRANSITION 

As has been discussed above, we have not seen 
the effects of small eigenvalues for P > Pc in our 
ma — 0.01, full QCD simulations to the extent 
expected, for example, from a model of dilute in- 
stantons. In an attempt to enhance such possible 
effects, we have repeated our calculation of (CC) 
for the pure gauge theory at and above the region 
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Figure 6. The quantity ((^C) plotted versus ■m(^ for 
P = 5.327 and ma = 0.025, It is compared with 
ma = 0.01 results for three values of /3. There 
is close agreement with the ma = 0.01, 13 = 5.3 
curve. 



of the dcconfining phase transition. The results 
are summarized in Figure |^ where we plot the 
"quenched" chiral condensate, {^Q as a function 
of j3 for a fixed mass value, = 0.001. 

The results are quite surprising. Because of the 
lack of symmetry for the Dirac operator, {QC,) 
shows dramatically different behavior in the three 
phases for f3 > Pc- The two phases where the 
Wilson line has a complex expectation value (here 
called complex phases) show the same behavior 
for {(Q because they are related by complex con- 
jugation. This complex phase shows spontaneous 
breaking of chiral symmetry, with values of (^C) 
that appear nearly continuous across the decon- 
finement phase transition. In contrast, the real 
phase has a vanishing value of (CC) in the limit 
TTiij —>■ and looks very much like the chirally 
symmetric phase seen in our full QCD simula- 
tions. 
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Figure 7. The quantity {^Q evaluated at 
= 0.001 is shown as a function of (3 for pure 
SU{3) gauge theory. The first-order, deconfining 
transition for Nt — 4 occurs at f3c = 5.6925. The 
multiple values for {(Q shown for [3 > Pc corre- 
spond to the real and complex phases with the 
smaller value of (CC) occurring in the real phase. 



4.1. Real phase 

In Figure ^ we show {(,(,) in the real phase as 
a function of quark mass, m^. One sees a frac- 
tional power dependence (CC) — 0.99(3) m° '''^'^\ 
very much like the behavior seen for /? > /3c in the 
full QCD simulations. As can be seen by compar- 
ing the results from the 16^ and 32'^ volumes, this 
behavior appears to characterize the infinite vol- 
ume limit. 

Figure ^ compares the dependence of {QC,) 
in this real phase with that of our full Nf = 2, 
ma=0.01 simulations. The results are somewhat 
similar, showing perhaps a 15% variation. How- 
ever, in contrast to the 0.01-0.025 comparison dis- 
cussed earlier, there is a clear systematic differ- 
ence between the quenched and full QCD results. 
As might be expected, the quenched results show 
an enhancement for small eigenvalues not seen in 
the ma ~ 0.01, two-flavor calculation. Never- 
the-less, we can describe the major effect of the 
fermion determinant on QCD thermodynamics as 
selecting the real phase (the phase with fewer 
small eigenvalues) as the physical phase. 
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Figure 8. A comparison of {((,) as a function of 
in the real and complex phases for f3 = 5.71 
and 16'^ and 32^ volumes. The real phase shows 
no perceptible volume dependence while the com- 
plex phase changes significantly. 



- 5.71 le^S quenched 

- 5.285 ma=0.01 

- 5.3 ma=0.01 

- 5.34 ma=0.01 



I 



5 * 



I 

! 



4.2. Complex phase 

We see chiral symmetry breaking in the com- 
plex, deconfined phase as might be expected from 
a dilute instanton gas model. In Figure ^ we ex- 
plicitly compare the dependence of ((fC) on quark 
mass for 16^ and 32^ volumes. As can be seen, 
the presumed finite- volume "knee", where (CC) 
abruptly begins to decrease with decreasing quark 
mass dramatically shifts to smaller quark mass on 
the larger volume. The shift by roughly the factor 
of 8 by which the volume has increased is exactly 
what is expected for such spontaneous symmetry 
breaking. 

The variation of ((^^) with is quite similar to 
that seen for /3 < /3c in our two-fiavor simulations. 
For smaller /3 a constant plus linear fit to (CC) as a 
function of works quite well but as /3 increases 
and the constant intercept decreases, a constant 
plus fraction power fit works significantly better. 
This effect is best seen on 32'^ volumes and in Fig- 
ure 10 we show (CC) versus for a series of /3 
values. This similarity with our (3 < f3c, Nf = 2 
simulations suggests the possibility of a second, 
chiral symmetry restoring transition in (CC) for 
this complex, deconfined phase. In Figure [l^ we 
plot the TO^ values from the constant plus 



Figure 9. A comparison of (CC) plotted as a func- 
tion of TTif computed the real, quenched phase 
at P = 5.71 (open circles) with our ma = 0.01, 
two-flavor results for a series of values of f3. Al- 
though the /3 = 5.34, full QCD results (solid in- 
verted triangles) match most closely for larger 
quark masses, as to^ decreases, the full QCD chi- 
ral condensate falls more rapidly. 



fractional power fits, versus (3. Clearly we cannot 
rule out such a second phase transition although 
a continuous, rapid decrease in (CC)|mf=o as f3 in- 
creases above the first-order critical point appears 
the most natural interpretation of the data. 



4.3. Inconsistent N 



f 



1 Phase Transition 



Combining our Nf — 2 results with this Nf = 
calculation suggests that the one-flavor phase 
transition may be quite different from that ex- 
pected on theoretical grounds, at least for cor- 
respondingly coarse lattice spacings. Since the 
density of small eigenvalues found in the "real" 
deconfined phase for pure QCD is insufficient to 
support a non-zero value of (CC)i the Nf = 1 
high temperature phase, with increased suppres- 
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Figure 10. The quenched chiral condensate in the 
complex phase, plotted versus using only 32'^ 
volumes for a series of values oi (3 > (3c- The 
curves are well fit by a constant plus a fractional 
power for m<; > 10"'*. 



sion of small eigenvalues coming from the fermion 
determinant, should also show exact chiral sym- 
metry with {xx) = in the limit m —>■ 0. Thus, 
for Nf = 1, {xx)m=o is non-zero for strong cou- 
pling but may well vanish for weak coupling which 
would require that Nf — 1 will show a phase tran- 
sition. This is in conflict with the usual analysis 
which concludes that there is no transition be- 
cause of the anomalous breaking of the theory's 
C/(l) chiral symmetry. 

The failure to see the behavior expected for 
Nf = 1 would be a serious cause of concern either 
about our understanding of continuum thermody- 
namics or about the coarse lattice spacings used 
in present-day lattice calculations. 

5. CONCLUSION 

It should be noted that we have adopted a new 
approach to the study of the axial anomaly using 
staggered fermions. Direct investigation of the 
axial anomaly is normally impeded by the ab- 
sence of a conserved, flavor-singlet axial current 
defined on the lattice. For example, a quantity 
like the rj' mass will receive contributions from 
both chirally-asymmetric lattice artifacts as well 
as the small Dirac eigenvalues of physical interest. 
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Figure 11. The constant term in our constant plus 
fractional power fit, plotted against /3. We cannot 
distinguish with certainty between the presence of 
a second, chiral symmetry restoring transition for 
(3 « 5.75 and a rapid decrease of (CC)l"^c ^ as 
P rises, although the latter interpretation appears 
to be favored. 



By using operators with a different number of fla- 
vors than appear in the Dirac determinant (oper- 
ators which obey an exact chiral symmetry on the 
lattice), we can study directly the effects of the 
relevant Dirac spectrum — a quantity determined 
solely by the number of flavors in the fermion de- 
terminant. The zero- momentum two-point func- 
tion u! of Eq. examined for Nf — 2 and the four- 
flavor chiral condensate (CC) discussed for Nf = 1 
are examples of this strategy. 
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